Dynamics of Two-Component Bose-Einstein Condensates Coupled with Environment 
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We investigate the dynamics of an open Bose-Einstein condensate system consisting of two hy- 
perfine states of the same atomic species which are coupled by tunable Raman laser. It is already 
suggested that the detuning between the laser frequency and transition frequency affect significantly 
on the dynamics of the pure condensate. Here we show that the detuning effect is suppressed by 
noise and dissipation caused by the environment. The increase of coherence and purity are also 
displayed for specific parameters. As a verification to the lowest-order approximation we derive the 
hierarchy of motion equations in the second-order approximation. It turns out that the former one 
can describe the dynamical evolution qualitatively for weak noise and dissipation and quantitatively 
for strong noise and dissipation. 

PACS numbers: 03.75.Lm, 03.75.Mn, 03.75.Gg, 03.65.Yz 



I. INTRODUCTION 

With the rapid experimental progress in manipulation 
of quantum gases multi-component atomic gases have 
now been an active research field in cold atomic physics. 
Since Myatt et al. first produced a binary mixture of 
condensates consisting of two hyperfine states {\F = 
l,m/ = -1) and \F = 2,m/ = 2)) of **^Rb fl|, the two- 
component atomic bosons of *^Rb in |-F = — ^1) 
and \F = 2^mf = 1) states received intensive study in 
experiments subsequently ■ variety of dynamical 
behaviors have been observed both in the Bose-Einstein 
condensate (BEG) and in the noncondensed sam- 
ple [5] . Meanwhile, theoretical works have predicted that 
the two-component Bose gas rnay exhibit exotic ground 
state and vortex structures |8l4l2l| as well as interesting 
dynamical properties, such as the Rabi oscillation of pop- 
ulation between the two states nonlinear population 
dynamics jT^ and quantum self-trapping [l4| . 

In the two-component atomic gas, the two hyperfine 
states (|F = = —1) and = 2,to/ = 1)) can be 

coupled by a two-photon transition which converts the 
^''Rb atom from one state to the other 0-11, Q. This 
transition thus drives the population dynamics. In this 
case, the two-component system is analogous to the Bose 
gas trapped in the double-well and both of them have 
also been investigated intensively in recent years [H- 
[Tg} . The dynamical properties of the condensate sat- 
isfy the time-dependent Gross-Pitaevskii (GP) theory be- 
cause the Bose atoms condensed in the ground state and 
the system can be described by the macroscopic wave- 
function. 

However, under experimental conditions, the conden- 
sates should be regarded as an open system coupled with 
the environment since the condensate usually coexists 
with noncondensed thermal cloud and atom loss is also 
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unavoidable. In this situation, noise and dissipation may 
play important roles. Many theoretical works have dis- 
cussed the noise and dissipation-induced effects in the 
open double- well condensates recently [13413 j ^'^'^ ii^" 
dicated that negligible changes in dynamical properties 
of the condensate appear, such as the decay of quan- 
tum self-trapping [l3|, decoherence [3| and dephasing 
[20j . The dephasingphenomenon has even been ob- 
served in experiments *2l||. However, sometimes dissipa- 
tions may result in enhancement of the quantum effect. 
For example, a thermally enhanced quantum-oscillation 
was observed during the domain formation process in the 
ferromagnetic spin-1 condensate [l^. Li et al. also found 
that the maximum spin squeezing can be reached even in 
presence of particle losses [23|. Most recently, Witthaut 
et al. showed that the dissipation could lead to enhance- 
ment of coherence in the double- well BEGs under specific 
conditions [23 |. 

In this paper we will investigates the open two- 
component BEGs. Although this system can be mapped 
into a double-well condensate model, it is of interest in 
its own right. For example, dynamics of this system can 
be controlled by tuning the laser frequency by which the 
two internal states of atoms are coupled O, [13] • So 
our main purpose is to discuss the noise and dissipation- 
induced effects when the two components are coupled 
by the detuned laser. The Hamiltonian is expressed us- 
ing the pseudo-angular-momentum operators and the dy- 
namical properties is evaluated by a hierarchy of ordinary 
differential equations of Bloch vectors, which are formu- 
lated on the basis of master equation method [3] ■ 

The paper is organized as follows. Section II formu- 
lates equations of motion in the first-order approximation 
for an open two-component Bose condensate with noise 
and dissipation. Section HI discusses dynamics of Bloch 
vectors for both the closed and open systems, respec- 
tively. Section IV is devoted to evolutions of the coher- 
ence and purity for certain given parameters. In Section 
V, we extend the motion equations to the second-order 
approximation and compare the results with those in the 
first order approximation. A summary is given in the last 
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section. 



II. MODEL AND METHOD 

The two-component BEC is described by the Hamiho- 
nian of second quantization in the single mode approxi- 
mation (natural unit ?i = 1 is used throughout the paper) 



H — ^ Hj + Hint + Hf , 
i=i,2 



(1) 



where 



with ujj = J d^rcj)* (r) -^^V^ + Vj (r) (f>j (r). The in- 
teraction term between the two hyperfine states of atoms 
takes the form of 



Hi, 



and the transition term is 



Hf^G 



Here gj (j=l,2) and gi2 are the effective intra- and 
inter-components interaction constants that can be con- 
trolled experimentally by tuning the atomic s-wave scat- 
tering length with Feshbach resonance techniques. The 
transition probability between both components with a 
small detuning is denoted by the Josephson coupling 
strength G. In the rotating wave approximation the 
phase Lp (t) ~ At, where A is the detuning between cou- 
pling laser frequencies and transition frequency \uji — uj2\ 
between two components. 

The Hamiltonian can be reexpressed by introducing 
the Schwinger pseudo-angular-momentum operators de- 
fined as = i {^a\a2 -f 0.201^, Ly ~ ^ (^a\a2 ~ aj^i^, 

Lz — ^ (^a\di — 0,20,2^ with Casimir invariant = 

T (t -'^)' '^here N — a\ai + a\a2 is the total bo- 
son number operator. L± = i: iLy and Lz sat- 
isfy the usual angular momentum commutation relation: 



Lz,L± = ±L-|-, and L-\-,L^ 
Hamiltonian can be written as 



2Lz. Therefore the 



H — ujqL. 



+ qLl + G(L+e"^^'^ +L^e-"^^'^^ (2) 



with Wo = uJi — UJ2 + {N — 1) {gi — 32) /2 and q — 
{gi + 92) /2 — gi2- In the Bose gas of \F — 1,to/ = —1) 
and |i^ = 2,m/ = l) states of ^''Rb, the effective in- 
teraction constants are known to be in the proportion 
91 -912 : ^2 = 1-03 : 1 : 0.97 [1]. These parameters can 
be tuned by magnetic Feshbach resonance [1^ so that 
they can be more different. 



With the Schwinger representation the general solution 
of the Schrodinger equation ih ^^^J^^'^ = H\tp{t)) govern- 
ing the dynamics of system can be obtained with a time- 
dependent unitary transformation and when the sys- 
tem couples with the environment, its dynamics will re- 
serve to the master equation 



H,p 



(3) 



where 
2 Z^ i 



p is the density operator, £p = 
2 2 ("-jP + /"^j ~ 2rijpnj) denotes the noise 

term and La = -\ Z]j=i,2 la-j + - Ihjpa}^ 

is the dissipation term. Accordingly, 7^ and 70 . 
describe, respectively, the strength of noise and the 
dissipation velocity for the component j. The one-time 
average of a system operator A can be calculated by 

i\ = tr 



Ap (t) . Here we define the single-particle 

Bloch vector and particle number as Sj (t) = 2tr Ljp (t) 
and n{t) = tr [(fii -|- /5 (t)], while the time deriva- 
tives of them are defined as Sj (t) = 2tr Ljp{t) and 

h (t) = tr np (t) 

Insert the Hamiltonian H into Eq. ([3|), we will 
find that the first-order operators Lj depend not only 
on themselves, but also on the second-order moments 

LiLj^. Similarly, the time evolution of the second-order 

moments depends on third-order moments, and so on. In 
order to obtain a closed set of equations of motion, the hi- 
erarchy of equations must be truncated at some stage by 
approximating the A'^th order expectation value in terms 
of all lower-order moments. Here we take the first-order 
approximations, i.e., l^LiLj^ k. (j^i^ {^■')' ^'^'^ thus the 
equations of motion for the Bloch vector take the formu- 
lation of 



Sx = -uioSy ~ qSySz - 2GszSmip (t) - Sx 
Sy = ujqSx + qsxs z - 2Gsz cos (t) - T^^Sy, 
Sz = 2G {sx sin ip (t) + Sy cos ip {t)) — T^^Sz - 



(4) 



Here the damping parameters Ti, T2 and the relative 
dissipation velocity fa are expressed as 



7a2 + 7ai 



The dissipation can be controlled artificially by shining 
a laser beam onto the condensates or be achieved by a 
forced radio frequency transition to an untrapped mag- 
netic substate 26]. The relative dissipation rate fa can 
be changed from to 1 and we take a medium value 0.5 
in the present paper if it is not given particularly when 
there exists atom loss. 

In this work we will investigate the phase coherence 
between two species and purity of the system, both 
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FIG. 1: The evolving Bloch vector for the closed system of 
iV = 1000 and A = 0. (a) g=0.001, a;o=0.0; (b) g=0.005, 
a;o=0.0; (c) g=0.008, cjq^O.O; (d) g=0.001, cjo=1.0. 



FIG. 2: The evolving Bloch vector for the closed system of 
= 1000 and ct;o=0.0. The initial state is chosen as Sx{0) = 
0.67V. (a) g=0.0, A =0.05; (b) g=0.0, A=0.1; (c) g=0.001, 
A=0.05; (d) q =0.01, A=0.05. 



of which are related with the Bloch vector by a (t) — 

2 (^d\a2^ I {fii +fi2) = y/s^ + s^/n and p(t) = |sp/n^ 

such that the dynamics of them can be investigated 
by solving Eq. Q for the definite initial condition. 
For convenience in the following evaluation the initial 
conditions s(0) = {sx{0), Sy{0), Sz{0)) will be taken as 
N{a, 0, Vl — o?) with N being the initial particle num- 
ber in the system. Sx{Q) = and Sz{0) = N {a—0) cor- 
responds to the case that all atoms occupy in the same 
hyperfine state while Sx{0) = N and s^(0) = (a=l) 
corresponds to the case that atoms populate equally in 
two hyperfine states. In the following evaluation we will 
consider the system of A'' = 1000. 



III. DYNAMICS OF BLOCH VECTOR 

In the present Schwinger representation all interesting 
physical quantities can be formulated with Bloch vec- 
tor. For example, for particle number conserved system 
coherence is the module of the projection of Bloch vec- 
tor on the Sx-Sy plane, its purity is the module square 
of Bloch vector and Sz{t) corresponds to the population 
imbalance between two components. So the dynamical 
evolution of Bloch vector is worth to study very much. 
This section will focus on the evolving Bloch vector and 
the explicit investigation of coherence and purity will be 
given in the next section. 



A. The Closed System 

Firstly, we investigate the dynamics of the system 
without coupling with environment, i.e., there is not noise 
and particle dissipation. In Fig. 1 (a-c) we display the 



Bloch vector for different nonlinear constants q for the 
case of detuning A = 0. The Bloch vector evolves peri- 
odically in the surface of Bloch sphere and forms a closed 
orbital, while different initial condition corresponds to re- 
spective orbital. 

For the weak nonlinear constants the orbitals behave 
as an almost perfect circular orbital parallel to the plane 
of Sx = and the fixed point belong to the line {sy — 
Sz — 0). That is to say, atoms shall oscillate periodi- 
cally between these two hyperfine states and in each pe- 
riod the average occupation probability is equal in both 
states. With the increase of nonlinear term the original 
circle orbital shall deform and their centers moved from 
the position oi Sz — {q = 0) to the north pole oi Sz = I- 
The fixed points move to the upper sphere. For example, 
Sz{t) approaches to 1.0 for all time as g = 0.008. This is 
referred to macroscopic quantum self-trapping (MQST) 
that atoms will prefer to stay in the same hyperfine state 
instead that during an oscillation period atoms occupy in 
equal probability in both states. The transition between 
periodical oscillation and MQST corresponds to the or- 
bital similar to lemniscates spreading in the surface of 
Bloch sphere. 

It deserves to notice that in Ref. Q those new formed 
two trajectories and fixed points F± for the strong in- 
teraction correspond to izo, respectively, while in the 
present evaluation we only choose the positive zq. In 
addition the transition from the periodical oscillation to 
MQST can be induced not only by the increased inter- 
action but also the increased initial imbalances between 
two states. 

The effect of ujq can be displayed by comparing Fig. la 
(a;o=0) and Fig. Id (a;o=1.0). It is shown that its value 
shall change the orientation of the orbital plane. The 
increase of ujq leads to the turn to the north pole of the 
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FIG. 3: The evolving Bloch vector for the system of A*' = 
1000 subject to noise for cjq^O.O, 7p=1.0 and Tf ^=0.0. (a) 
q =0.00, A=0.0; (b) g=0.005, A=0.0; (c) g=0.001, A =0.05; 
(d) g=0.005, A=0.05. 



FIG. 4: The evolving Bloch vector for the system of A'' = 1000 
subject to dissipation for cjo=0.0, 7p=0.0 A=0.05 and T^^ 
= 1.0. (a) g=0.001; (b) g=0.005; (c) g=0.008; (d) g=0.01. 



orientation of orbital plane, which is also an expression of 
MQST. In a word, for the closed system of nondetuning 
MQST can exhibit by tuning 9 ^^nd initial imbalance, 
whose orbital center deviate from ~ 0. 

In Fig. 2 we display the orbitals for the situation where 
there exist detune between the transition frequency and 
the coupling laser. In this case the Bloch vector will no 
longer evolve periodically and its path behave as a spiral 
line in the surface of Bloch sphere. With the time delay 
the path of Bloch vector wind each other and looks like a 
basket when the evolving time is long enough. By com- 
paring Fig. 2a with Fig. 2b and Fig. 2c with Fig. 2d we 
found that the sparser the net grid, the bigger the detun- 
ing and the denser the net grid, the stronger the nonlinear 
constant. The denser track means that the evolution of 
Bloch vector and thus the dynamics of two-component 
BECs is more close to the pure periodicity in short time. 
In spite of the complicated orbital its projection on Sz 
axis is still periodical oscillation while the period is larger 
a little than the case without detune. The evaluation 
show that its period become longer although the atoms 
still transit periodically between two states and the tran- 
sition point between delocalization and self-trapping also 
deviate slightly from the case without detune. 

B. The System Coupling With Environment 

After investigating the dynamics of closed system, we 
focus on the system coupling with environment. The 
effects of noise and dissipation are shown in Fig. 3 and 
in Fig. 4, respectively. Here different orbitals in Bloch 
space correspond to different initial condition. 

For the system subject to noise the module of Bloch 
vector will not be conserved rather decrease to zero in 



a spiral during the dynamical evolution. It would not 
be on the surface of Bloch sphere rather shrink to the 
original point of coordinate frame even for very weak 
noise, whose strength only affect the evolving time to 
the original point. That is to say, when suffering noise 
instead of dissipation the atom number populated in two 
components would always tend to balance. In addition 
the purity of system tend to zero that corresponds to the 
module of Bloch vector. This will be discussed further 
in the later section. By comparing Fig. 3a with Fig. 3b 
and Fig. 3c with Fig. 3d we found that the stronger 
nonlinear interaction constant and large detuning shall 
induce that the Bloch vector run as an irregular orbital. 

In Fig. 4 we plot the orbitals for the case that there are 
only dissipation instead of noise for different nonlinear in- 
teraction constants q. The Bloch vector shall be always 
on the surface of Bloch sphere but the orbitals become 
irregular and not closed. In this situation for different ini- 
tial atom population in two components the final state 
shall be indefinite completely. With the increase of dissi- 
pation the atomic loss rate become faster and the lifetime 
of BECs become shorter and shorter. This means that 
along with the BECs oscillation between two hyperfine 
states aperiodically, more and more atoms escape from 
the system. As the interaction constant increases the 
shape of orbitals become more and more complicated. 

Fig. 5 displays the evolution of Bloch vector for the 
system of 1000 atoms and q — 0.005 subject to both 
noise and dissipation. Here /a=0.5 in Fig. 5 (a)-(e) and 
/a=0.0 in Fig. 5f. It is shown that for different initial 
states it will always evolve to the same point in the Bloch 
space, which correspond to the same final state, while the 
change of noise and dissipation shall determine its evolv- 
ing path and dissipation velocity to the final state. In 
addition the joint effect on the BECs induces that the 
evolving path to the final state become more simple than 



FIG. 5: The evolving Bloch vector for the system of TV = 1000 
subject to both noise and dissipation for a;o=0.0 and g=0.005. 
(a)-(c) 7p=2.0, rf^ = 1.0; A=0.0 (a), 0.05 (b), 0.1 (c). (d) 
7p=1.0, Tj-^ =2.0, A=0.1. (e) 7p=2.0, T^^^2.0, A=0.1. (f) 
7p=2.0, r-'=2.0, A=0.1, fa^O.O. 



the system subject to only noise or dissipation. Accord- 
ing to Fig. 5 (a)-(c) (A=0.0, 0.05, 0.1) the detuning 
between transition frequency and the frequency of laser 
induce no explicit effect. According to Fig.Se and Fig. 
5f, the relative dissipation rate of each components shall 
result in different final position in Bloch space, i.e., dif- 
ferent final states. 

To sum up, noise always induce that atoms populate in 
two components evenly and atom dissipation determine 
the lifetime of condensate. Both noise and dissipation 
shall result in the irregular evolving path of Bloch vector, 
i.e., the aperiodic dynamics. In this situation the system 
always arrive at the same final state and the evolving 
path shall be more simple compared with the case subject 
to only noise or dissipation. 



FIG. 6: The coherence dynamics for the system of q=0.001, 
a;o=0, Go=l, A=0 and /„ =0.5. (a) 7=0.0, rj-^=0.0; (b) 
7=1.0, T-^=0; (c) 7 =0.0, r-^=2.0; (d) 7=2.0, Tj"^=3.0. 



1.0 
0.8 
0.6 
0.4 
0.2 






s_(0)=0.0 




s^(0)=0.5 




s_(0)^0.9 




s.[0)=1.0 ■ 




FIG. 7: The time evolution of purity for the system of 
q=0.001, a;o=0, Go = l, A=0 and /„ =0.5. (a) 7=0.0, 
rr'=2.0; (b) 7=2.0, r-i=0; (c) 7 =2.0, (d) 7=2.0, 

T^^=3.0. 



IV. COHERENCE AND PURITY 

In this section we will study the evolution of coherence 
and purity for different noise and dissipation. 

In Fig. 6 the evolution of coherence is displayed for 
different initial conditions. It is shown that for a closed 
system the coherence will always oscillate periodically as 
long as the initial population difference in each compo- 
nent deviates from zero (a 7^ 1) while it will preserve 
the strongest coherence when the atoms initially popu- 
late equally at two components (a — 1) (Fig. 6a). When 
the BECs suffer noise instead of dissipation the coher- 
ence will decrease oscillationly (Fig. 6b) , and when there 
exist dissipation instead of noise the coherence still os- 
cillates quasi-periodically (Fig. 6c). For the system ex- 
posed to both noise and dissipation the coherence show 
rich evolving dynamics determined by the initial popula- 
tions (Fig. 6d). For the case of equal population initially 



(a=l) the coherence between two components decrease 
firstly and then increase, which is similar to the stochas- 
tic resonance. For the polarized case initially (a=0) the 
coherence show a process of first increase from zero to 
the strongest value and then following a stochastic res- 
onance. While for the initial population between these 
two limits (0 < a < 1) coherence will show the evolving 
process of first increase or decrease and then following a 
stochastic resonance. 

In short, particle loss do not temper with the coher- 
ence dynamics greatly and the obvious effect of noise is 
the reduction of the maximum coherence. So by tuning 
the strength of noise and dissipation rate we always can 
observe the process of coherence enhancement. The sim- 
ilar effect on the purity can also be seen (Fig. 7), where 
dissipation show no effect on purity. 

The evolving purity is shown in Fig. 7. According to 
Fig. 7a (7 = 0) and Fig. 7b (Tf ^ = 0) the dissipa- 
tion shall not bring about the change of purity as long 
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as the noise term is equal to zero and the noise shah 
result in the decrease of purity. By tuning the relative 
strength of noise and dissipation rate the purity increase 
again after rapid decrease. It turns out that for suitable 
parameter the purity can be enhanced to a very large 
value. For instance, it reach 0.7 for 7=2.0, r|~^=3.0. 
In addition, the purity evolving do not show qualitative 
difference for different initial condition. It deserves to 
notice that the relative dissipation rate fa play a criti- 
cal role on whether the enhancement of coherence and 
purity can be displayed. Generally, the bigger fa shall 
result in the larger coherence and purity so we take a 
medium value of /q = 0.5 here. 




FIG. 8: Comparison between the lowest and second order 
mean-field approximation. q=0.001, cJo=0, Go = l, A=0, fa 
=0.5. (a) 7=2.0, ri"^=3.0; (b) 7=0.5, rj"^=0.5. 



V. SECOND-ORDER APPROXIMATION 

In order to check the accuracy of the lowest-order ap- 
proximation, we now take into account the second-order 
correction and compare the results from both cases. Up 
to the second-order approximation, the hierarchy equa- 
tions are truncated by approximating the third-order ex- 
pectation value < LiLjLk > as the following: 



LiLjLk 



LjLk 



L. 



2 L, ) L 



Thus the hierarchy of motion equation shall be composed 
of the Bloch vector and the second-order moments An — 



4(< LiLj + LjLi > 
formulated as below 



-2 < Li >< Lj >). They can be 



-UJQSy 



2sySz) — 2Gsz sin ip (t) 



(A^2 -t- 2sxSz) - 2Gsz cos (p (t) 



T2 Sy 

-1 



2G sin ip {t) + Sy cos ip (t)) - - faU, 

4Gsmip{t) - 2q [A^zSy + A^yS^] 



-2woA 



xy 



-2T-^Axx + 2 {T-^ - T-^) {Ayy + 2^) + 2T-\, 
2ujoAj:y - AGcosip{t) Ayz + 2q [A^ySz + Ay^Sx) 



-2T^^A 



vv 



2 {T,-' - Tf 1 



A, 



2Tf 



Azz 



Axz = 



A„ 



4G (cos ip (t) A. 

1^0 {AxX — Ayy) - 
+ q{AxxSz - Ayy 

-2 {2T-'^ - T-'^' 
-UJoAyz - qAyzS 
-2Gsuvip{t) [Az. 
loqAxz + q {AxzS 
~Az 



^ 2s 

-sin(^(t) Axz) ~ 2T-^ [Azz - n] 
2G cos ip (t) Axz ~ 2G sin ip (t) Ayz 

z Axz^x A 



Sj,) 



Ax,. ~ 4 {T-^ ~ T^^) 



zz'^y 



2G COS If (t) Axy 



xy 

qAz 

Axx]- [T2'+Ti^] ^xz 

AzzSx) + 2G COS ip{t) (A, 
A„ 



- 2Gsin^ (i) Axy - [T^^ + Tf j ^yz. 

Fig. 8a and 8b display the coherence a and purity 
P obtained based on both the lowest-order and second- 
order approximation. For strong noise and dissipation 



case, as shown in Fig. 8a, there is not obvious differ- 
ence between the lowest-order and second-order results. 
For weak noise and dissipation case, although the sec- 
ond correction becomes more and more significant with 
the time, the tendencies of the curves are similar [see 
Fig. 8b]. Therefore, we conclude that the lowest order 
approximation are creditable qualitatively in the full pa- 
rameter regime, especially in the presence of the noise 
and dissipation. 



VI. SUMMARY 

In conclusion, we have investigated dynamics of the 
open two-component EEC in the Bloch representation. 
The two components correspond to two hyperfine states 
of the same specie of atom which are coupled by the 
detuned Raman laser. We have calculated the evolving 
Bloch vector, coherence and purity for both the closed 
system and open system. 

For the closed system, when transition frequency 
match with Raman laser, the evolving Bloch vector form 
a closed orbital in the surface of Bloch sphere. This 
mean that all physical quantities for this system evolve 
periodically. As the nonlinear term becomes strong the 
MQST exhibits. In the detuning case, the orbital shall 
not be closed and wind into a basketlike path. In com- 
parison with the non-detuning case, although atoms still 
transit between two states periodically, its transition pe- 
riod become longer. For the system subject to noise the 
atoms always evolve to a balance state that corresponds 
to the zero point in the Bloch space. In the situation of 
strong nonlinear term and large detuning the orbital of 
Bloch vector shall be irregular. For the system of dissi- 
pation, along with atoms oscillate aperiodically between 
two states more and more atoms escape from the conden- 
sates. The evolving Bloch vector behave as an irregular 
orbital in the surface of Bloch sphere. When the conden- 
sates subject to both noise and dissipation the system al- 
ways evolve to the same final state although the evolving 
path might be different for different noise strength and 
dissipation rate. It deserves to notice that the detuning 
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between transition frequency and laser frequency change 
obviously the dynamics of the closed system, while the 
detuning effect is not displayed obviously in the presence 
of noise and dissipation. 

We also studied the coherence and purity for both 
closed and open system. For the closed system the co- 
herence always oscillates periodically. As subject to noise 
the coherence shall decay to zero during the period of evo- 
lution, while the particle loss induce no obvious effect on 
the coherence. By tuning the noise strength and dissipa- 
tion rate the interesting effect of coherence enhancement 
exhibits. The investigation of purity show that the sys- 
tem preserve the maximum purity as long as it is not 
subject to noise even if the dissipation exists. In spite 
that the noise shall decrease the purity, it can be en- 
hanced to a large value by tuning the noise strength and 
dissipation. 



Wc discussed the second-order correction of the mas- 
ter equations. Comparison between the results from the 
lowest-order approximation with those from the second- 
order one has shown that the former can describe quali- 
tatively the dynamics of both open and closed BECs and 
with the enhancement of coupling with environment its 
quantitative result shall be exact. 
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